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The ac Stark shift in the presence of spontaneous decay is typically considered to induce an 
effective dephasing with a scattering rate equal to r s |£J| 2 /A 2 , where T s is the spontaneous decay 
rate, Ul is the laser transition coupling, and A is the detuning. We show that under realistic 
circumstances this dephasing rate may be strongly modifed due to non-Markovian behavior. The 
non-Markovian behavior arises due to an effective modification of the light-atom coupling in the 
presence of the ac Stark shift laser. An analytical formula for the non-Markovian ac Stark shift 
induced dephasing is derived. We obtain that for narrow laser linewidths the effective dephasing 
rate is suppressed by a factor of Q 2 , where Q is the quality factor of the laser. 

PACS numbers: 03.75.Dg, 37.25.+k, 03.75.Mn 


Introduction — The ac Stark effect, or light shift, is 
one of the fundamental phenomena in atom-light inter¬ 
actions El- The basic effect results from applying an 
off-resonant laser between two atomic levels, which cre¬ 
ates an energy shift equal to |S1| 2 /A, where Ul is the laser 
transition coupling and A is the energy detuning [2|]. The 
effect is responsible for a wide variety of techniques in 
atomic, molecular, and optical physics. Since the energy 
shift is proportional to the intensity of the light, it is 
the basis of optical traps, where a spatially varying light 
field creates regions of potential minima [3|. Optical lat¬ 
tices created by two counterpropagating laser beams in 
a standing wave configuration creates a periodic poten¬ 
tial due to the same effect [dj. Laser cooling techniques 
such as Sisyphus cooling also rely on the ac Stark ef¬ 
fect by exchanging kinetic energy for potential energy 
in an optical lattice geometry |4j, |5j . Such trapping and 
cooling techniques are fundamental for quantum infor¬ 
mation applications such as quantum simulation [fi] and 
quantum computation jj| based on neutral atoms. The 
control of how much the light shift occurs is an impor¬ 
tant problem for atomic clocks, where magic wavelengths 
are used to ensure cancellations of the shifts [IHol • The 
observation of the ac Stark effect is now of fundamental 
importance to other quantum systems being investigated 
for quantum information applications, such as supercon¬ 
ducting qubits El jBose-Einstein condensates El, El, 
and quantum dots |14(. 

Under the Markovian approximation the ac Stark shift 
in the presence of spontaneous decay introduces a de¬ 
phasing (the “scattering rate”) on the ground states cou¬ 
pled to the laser according to El is 


r s |0 | 2 




A 2 ’ 


(1) 


where T s is the spontaneous decay rate. It is often as¬ 
sumed that under typical experimental situations that 



FIG. 1: (Color online) The ac Stark shift configuration. An 
atom is illuminated with an off-resonant laser between levels b 
and e. The laser is assumed to have a Lorentzian distribution 
with a central frequency ljo, maximum intensity |ao| 2 and 
linewidth A. The frequency dependent coupling between the 
levels b and e is g{cd), and the spontaneous emission rate is 
T s . The laser is off-resonant by a frequency A. a , b label the 
ground states of the atom. The coherence between the ground 
states is dephased due to the ac Stark shift of the laser. 


the Markovian approximation is sufficient to capture the 
basic effects of decoherence in quantum systems. How¬ 
ever, this assumption depends upon several factors, and 
in particular for atomic systems (the “system”), the ef¬ 
fects due to radiation field (the “bath”) are known to be 
either Markovian or non-Markovian depending on the cir¬ 
cumstances El, El- It has been shown that such effects 
are important to understand the dynamics as shown in 
several experimental systems EU, [l9j , and has been even 
shown to be controllable such as to enter both regimes 
[20j . Understanding the nature of non-Markovian dy¬ 
namics is naturally a very important topic for quantum 
information science, where the aim is to control aquan- 
tum system for use in technological applications [21H24I | . 
For this reason several works have studied the funda- 
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mental problem of a two-level system in a variety of non- 
Markovian settings [25l - [33j . Most works have focused 
on the modification of spontaneous decay in the non- 
Markovian regime. However, to our knowledge, no cal¬ 
culation of the non-Markovian scattering rate has been 
performed before. In this paper, we analytically calcu¬ 
late the ac Stark shift dephasing in the non-Markovian 
regime when illuminated by a narrow linewidth laser. 

The standard way to calculate the ac Stark shift scat¬ 
tering rate © is by assuming a master equation involving 
a ground and excited state, with the excited state hav¬ 
ing the possibility of spontaneous decay fl3 | (see Fig. [TJ. 
Specifically, the coherent part of the evolution is deter¬ 
mined by the Hamiltonian 

H/h = A|e)(e| + (H|6)(e| + H.c.) (2) 

in the presence of spontaneous decay described by a mas¬ 
ter equation 

l = < 3 > 

where the Lindblad decay term is C[0\p = 2 OpO^ — 
O^Op — pO^O. We show that this commonly assumed 
framework in fact misses several aspects which turn out 
to be rather important to the problem. The first is that 
the spontaneous decay rate T s is derived without the 
presence of the laser. The laser only enters the problem 
after the (0 is constructed. As we show, if the dephasing 
is derived in the presence of the laser, qualitatively differ¬ 
ent results are obtained. The second related issue is that 
when T s is derived in ©, the Markovian assumption is 
quite reasonable, as the coupling to the environment can 
be considered to be smooth. However, in the presence 
of a laser, the effective coupling no longer satisfies this, 
and non-Markovian effects become more important. We 
show in this paper that these aspects strongly modify the 
standard dephasing rate ©• In particular, under typical 
experimental parameters we expect that non-Markovian 
effects are in fact very important to consider. Remark¬ 
ably, for narrow linewidth lasers this results in dephasing 
rates that are greatly suppressed in comparison to the 
standard result © . 

Model Hamiltonian — Our derivation of the ac Stark 
shift induced dephasing proceeds as follows. We con¬ 
sider an excited state |e) and ground states of an atom 
| a,b), coupled to a continuum of electromagnetic modes 
according to 

H/h =w e |e)(e| + w a |a)(a| +ui b \b)(b\ + y] a > k pjpk 

k 

+ '^2{gk\e)(b\p k +H.c.) (4) 

k 

where p k , uj k are the respective annihilation operator and 
frequency of the fcth bath mode, io a ^,e are frequencies of 
the levels a, b, e, and gk describes the coupling to the elec¬ 
tromagnetic field via electric dipole approximation (see 


Fig. Q). We assume for simplicity that the photons only 
couple to level b, which may arise in practice due to selec¬ 
tion rules of the laser transition. Our results can be easily 
generalized as long as the bath modes for levels a, b may 
be treated independently. In the usual treatment one 
would eliminate the bath modes in © to obtain above 
the master equation under a Markovian approximation, 
without the presence of the laser. However, in our pro¬ 
cedure instead we will derive the dephasing master equa¬ 
tion in the presence of the laser. In the presence of the 
off-resonant laser, there will be a second order virtual 
process of the atom being excited to higher energy level 
and relaxing back to the ground state. For a large de¬ 
tuning A between the laser frequency u o and the atomic 
transition frequency, we can eliminate the excited state 
adiabatically. Setting the zero point energy at level b 
and following the standard adiabatic elimination proce¬ 
dure by setting i = H\e') = 0, where \e') = e luJot \e), 
we obtain within the rotating wave approximation 



where a z = |a)(a| - |6)(6|, w a b = w a -w b , A = u> 0 - (u> e - 
cob) and constant terms have been dropped. This Hamil¬ 
tonian describes the phase coupling occurring during the 
virtual transition to excited state and relaxing back to 
the ground state. 

The presence of the laser may be taken into account 
by displacing the photon operators with the Hamiltonian 
-ffiaser = J2 k F kPk+F k p[, where F k is the laser amplitude. 
The laser displacement may be eliminated by a transfor¬ 
mation q k = Pk + c^k = D~ 1 (ak)pkD(ak) 1 where a k is the 
resulting coherent state amplitude, D(a k ) = e akPk ~ akPk 
is the displacement operator [13] ■ The Hamiltonian is 
then 

H e{{ /h = + J2 UJ kqtqk-J2^~ crZ ^ akq k +R - c -) 

k k 

( 6 ) 

where u r = oj ab + J2k renormalized fre¬ 

quency between the ground states, and we have assumed 
that Wfc as the bath modes that are relevant to 

the atomic transition are much larger than the second 
order shifts in the energy. Moving to interaction pic¬ 
ture with respect to bath and the system, the interaction 
Hamiltonian can be written as 

Hi(t)/h %^‘( iUkt + H.c.). (7) 

k 

Non-Markovian scattering rate — The time evolu¬ 
tion of the bath and spin can be obtained by evolv¬ 
ing the total density matrix pr(0) with the time evo¬ 
lution operator U(t) from an initially unentangled state 
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p s (0) <S> ]~[fe \otk)k{a k \k in terms of the original variables 
(■ a,Pk )■ For the displaced basis (cr, <?*,), the initial state is 


PT(0)=P8(0)®Y[\0)k(0\k. ( 8 ) 

k 

In the interaction picture, the time evolution operator 
U (t) is given by 


U(t) = T^exp[-^J dsH^s)} (9) 

where T<_ is the time ordering operator. The above time 
ordered product is usually difficult to evaluate. However, 
the interaction Hamiltonian Hi in our case simplifies its 
evaluation. We note that [cr z (t),Hi(t')\ = 0 and the 
unequal time commutator 

[H^H^)] = -2 ( 10 ) 


which is constant in all operators as <f>(i) = 

sinwfc(i). This allows us to write the time 
ordered evolution operator in the form 


U(t) = exp 


x exp 


i dti dt 2 ^{t\ - t 2 )0(ti - t 2 ) 


a z ^ J^L(a fe f k (t)ql - H.c) 


( 11 ) 


where f k (t) = e * -1 . The first factor in U(t) is a global 
phase factor and does not involve any dynamical change 
to the time evolution of the system. 

We may now obtain the effective dynamics of the spin 
alone at a time t by taking the trace over bath degrees 
of freedom p s (t) = Tr b[U (t) pt (0)W (t)]. Taking p a s a ' = 
(<j|p s |cr') with u z \cr) = cr|cr), we obtain 

P7\t) = P r'(0)exp h (g ~ a,)2 r(t)] 

where 

T(t) = y- 4 lgfc| 4 |afc| 2 (1 ~cosw fc t) 

V 4 A 2 w 2 


( 12 ) 

(13) 


As expected the dephasing preserves the diagonal ele¬ 
ments of p s (t) and eventually reduces off-diagonal density 
matrix elements to zero. 

In most experimental situations the ac Stark shift is 
induced by a narrow linewidth laser, in which case the co¬ 
herent field distribution is well described by a Lorentzian 


ck 2 (w) 


Kl 2 A 2 

7T (w — Wq) 2 + A 2 


(14) 


Here A represents the spectral width, |«o| 2 is the inten¬ 
sity, and wo is the center frequency of the laser. The 
decoherence function (fl3l) may be evaluated by using the 
standard electric dipole approximation in free space to 


obtain an expression for gk as shown in the Appendix. 
The decoherence function then reduces to 


m = 


r s |fi | 2 


A 2 


A 2 t 


-Wg + A 2 


A[(wq — A 2 )(l — e xt coswot) — 2woAe xt sinwot] 


2(w 2 


A 2 ) 2 


(15) 

The coefficient of (fl5l) is precisely that given in the stan¬ 
dard expression for the scattering rate r« in Eq. (UJ. In 
terms of the dimensionless time r = Tm^, there are two 
non-trivial parameters in (1151) : the laser quality factor 
Q = wo/A, and the ratio R = A /Tm which allows us to 
rewrite the decoherence function as 

T(r) = 


Q 2 + l 


(' Q 2 — 1)(1 — e Ht cos (QRt)) — 2 Qe Kt s'ui^QRt) 


2 R{Q 2 +1) 2 


(16) 


Analysis of decoherence function — Let us first under¬ 
stand how (fl6l) reduces to the standard Markovian limit. 
Figure (2])a) shows the coherence factor e~ r ^ for various 
values of Q with R 3> 1. We see that the the Markovian 
curve is recovered for Q = 0 and R 1. Physically, this 
can be understood by looking at the role of the various 
parameters. The typical time scale over which there ex¬ 
ists an intrinsic time evolution of the system is ~ 1/wo- 
Meanwhile as A is the spectral width of the coupling, the 
bath correlation time ~ 1/A. The Markovian approxima¬ 
tion implies that system evolves over large time compared 
to very fast bath dynamics, implying Q <C 1 and R 1. 
In this limit (fTBl) reduces to V(t) = TMt, agreeing with 
standard Markovian dephasing 0. 

However, the Markovian limit is not always the physi¬ 
cally relevant regime. In practice the linewidth A of the 
laser is much narrower than the transition frequency wo, 
implying Q 1. Taking numbers for 8 7 Rb for exam¬ 
ple [33], the spontaneous emission rate is T s « 19 MHz, 
which may typically give Tm ~ MHz. The transition fre¬ 
quency is wo ~ 300 THz, and the linewidth in the range 
A ~ kHz to MHz, implying R < 1. These parameters 
suggest that the Markovian limit is in fact the opposite 
regime in terms of the magnitudes of Q and R to the 
physically relevant regime. 

The decay of the off-diagonal terms of the density ma¬ 
trix at various values of Q for the more physically rel¬ 
evant R < 1 regime is shown in Fig. 0b). It can be 
observed that as Q increases, there is a larger deviation 
from Markovian behavior. The general trend is similar 
for all values of R as can be seen on comparison to Fig. 
0a) . The term that is primarily responsible for this is the 
first term in (1161) . where the Q 2 + 1 term appears in the 
denominator. Considering that Q 1 is the physically 
relevant regime, this suggests defining another timescale 


r M r s |H| 2 A 2 

Q 2 A 2 Wq 


(17) 
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FIG. 2: (Color online) Decay of coherence with time at dif¬ 
ferent values of the parameter Q = uj o/A and R = A/Tm- 
Plots are for Q values as marked with (a) R = 100, (b) (c) 
R = 0.01, (d) R = 10 -5 . Plots (a) and (b) are in units of 

p q2 

the Markovian decay time Tm = ■ Dotted line shows 

the Markovian approximation e~ TMt . Plots (c) and (d) are 
in units of the ac Stark shift dephasing rate r ac = Vm/Q 2 ■ 
The dotted line shows the approximation e~ r ° ct . 


Figure [U(c) (d) shows the behavior of the coherence func¬ 
tion in terms of the rescaled time t' = Y ac t for two dif¬ 
ferent values of R. We see that for large Q the coherence 
function coincides with exponential decay in units of the 
rescaled time. 

Before reaching the limiting behavior of exponential 
decay with rate ED. we observe some oscillatory behav¬ 
ior Figure [2])c)(d). This occurs clearly due to the second 
term in m■ For Q > 1 we may approximate this as 


F(r') 


j (1 - e- flQV cos(11QV)) 
+ 2 RQ 2 


(18) 


In order for the oscillatory term to be visible, we require 
the exponential decay term in (fl8l) to be not too fast 
giving RQ 2 < 1, but simultaneously the oscillation fre¬ 
quency should be faster than the overall decay envelope 
RQ 3 > RQ 2 > 1. The strongest oscillations therefore oc¬ 
cur when RQ 2 ~ 1, which agree with the numerical plots 
in Figure 0(c) (d). For larger values of Q, the oscillatory 
term in (TT51) is suppressed, and we return to pure expo¬ 
nential decay with a rate m3. For typical experimental 
parameters as given above we have RQ 2 1, hence much 
lower values of Q are necessary to see the oscillatory be¬ 
havior. The deviation from an exponential decay can 
be attributed to the memory effects developed initially, 
typical of non-Markovian behavior. As RQ 2 = lUq/Tm^ 
it is natural to expect that such a parameter should be 
sufficiently large to observe oscillations, which may be 
achieved by having a narrow linewidth and low spon¬ 
taneous emission rate. In general, three time scales in 
an open system exist to characterize non-Markovian dy¬ 
namics: (i) the time scale of the system ~ 1/ujq (ii) the 


time scale of the bath ~ 1/A (iii) the mutual time scale 
arising from the coup ling between the system and the 
bath ~ 1 /Tm [30, .35, j36(. In this case the criterion for 
the strongest oscillatory behavior are satisfied when these 
are all approximately the same. 

Comparison to past works — Several works have stud¬ 
ied non-Markovian dynamics of a qubit with a Lorentzian 
bath before ]2Q, |27j, (3(|. Here we discuss whether the re¬ 
sults of these past works could be used to derive our 
primary result m- One potential approach to deriving 
this equation would be to use the standard expression 
for the ac Stark shift scattering rate ©> and replace the 
non-Markovian spontaneous decay in the place for T s . 
We find that such an approach does not yield the same 
results as our m- „ 

Taking the example of Ref. [27], there the non- 
Markovian decay of the density matrix p of a two-level 
system in a Lorentzian bath is derived. From Ref. [27} , 
the excited state density matrix element obeys 


Pee(t) = e 


—A t 


cosh(y<5) + isinh(y<5) 


Pee{ 0), (19) 


where <J = y 1 — and r s is the Markovian sponta¬ 
neous decay rate. In the weak coupling limit T s -C A, 
there is a very broad coupling to many frequency modes, 
which gives Markovian behavior. Here S « 1 — ^ and 
the decay function (flT)l) gives purely exponential behav¬ 
ior. To first order the decay function may be approx¬ 
imated as p ee (t) w e~ Tst l 2 p ee (Q), which is nothing but 
standard Markovian spontaneous decay. In this limit the 
ac Stark scattering rate would be exactly the same as the 
standard expression ©■ Meanwhile, in the strong cou¬ 
pling limit T s A, the linewidth of the laser is extremely 
narrow, which gives rise to strongly non-Markovian be¬ 
havior. Here we may approximate 5 = ~ and 


Pee(t) = e~ Xt 

which cor respond s to damped oscillations at frequency 
CInm = yAT7/2 and a decay envelope with rate A. Tak¬ 
ing r s = A and substituting into © we would obtain 

r ac = A ^J . This again does not yield the scaling factor 
Q 2 as derived in ED- 

The above approach and ours mainly differ in what 
stage the bath is eliminated to produce an effective mas¬ 
ter equation. It is helpful to examine what order effect 
we are examining in terms of the original couplings gt- 
In our derivation we start with an effective coupling of 
the form oc g^.(T z p\pk to the bath modes which directly 
gives the dephasing of the qubit. As is typical of all mas¬ 
ter equation derivations, the dephasing rate is the square 
of the coupling, hence our final expression is effectively 
an expression that is fourth order in r ac oc g^, as seen 
in m- On the other hand, if one uses the spontaneous 


cos (tl NM t) + \l 7 ^ sin (V,NMt) 


Pee (0) 
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emission as a starting point, the dephasing is not cal¬ 
culated directly, but is inferred from the non-Markovian 
spontaneous emission rate. In this case the bath is elimi¬ 
nated at first order in the coupling, resulting in the spon¬ 
taneous emission being a second order effect oc g\. 
Then if we substitute this into there is another fac¬ 
tor of g\ from the id 2 , hence both give expressions that 
are fourth order in g k . However, our derivation is a di¬ 
rect fourth order calculation, whereas © is combining 
two second order results. Combining two second order 
results does not yield generally the same results as a di¬ 
rect fourth order calculation, which is the case here as 
well. 

Conclusions — We have derived the ac Stark shift scat¬ 
tering rate in the fully non-Markovian regime. In con¬ 
trast to the standard procedure where the laser field is 
put in after the bath is eliminated, in our procedure the 
bath is eliminated only at the end to obtain the dephasing 
rate. While our results reduce to the standard expression 
for the scattering rate in the Markovian limit, we find 
that this is not the typical physically relevant regime. In 
the physically relevant regime, the scattering rate is bet¬ 
ter approximated by T ac (Eq. (TT71) 1 . which is equal to 
the standard Markovian scattering rate divided by the 
square of the laser quality factor Q , which can be quite 
large in practice. We note that in terms of the order of 
the result, both are fourth order in terms of the atomic 
transitions. Thus the difference does not result from a 
difference in order of the calculations, and is a result 
of a consistent non-Markovian treatment of the dephas¬ 
ing induced in the presence of the laser. Furthermore, 
CEB cannot be derived by simple methods using past 
works deriving non-Markovian spontaneous emission in 
a Lorentzian bath [2(1 [23, [3Cj • We attribute this to the 
fact that our result is a direct fourth order calculation, 
whereas © combines two second order effects. 

Physically, we may understand the difference as fol¬ 
lows. In the standard approach, dephasing may occur 
because during one of the transition of the ac Stark shift, 
instead of being driven back down by the laser, sponta¬ 
neous emission may occur. As the spontaneous emission 
in © is derived without the presence of the laser, the 
excited mode couples to the full range of frequencies of 
the bath modes. On the other hand, in our derivation, 
the presence of the laser linewidth bosonically enhances 
emission into the bath modes that are excited by the 
laser. Thus spontaneous emission occurs only at those 
modes that are excited by the laser, which is much nar¬ 
rower than that predicted by ©. For fewer bath modes, 
the dephasing rate is naturally much smaller, in this case 
by a factor Q 2 . 

In most practical applications of the ac Stark shift, 
the experimental parameters are chosen such that even 
assuming the rate ©, the detuning is made large enough 
that the scattering rate is made negligible compared to 
the light shift. The current result would imply that par¬ 


ticularly for narrow linewidth lasers, the detuning re¬ 
quirements could be made less, offering an alternative 
approach to controlling the dephasing. This would be 
important for applications where spontaneous emission 
is a serious drawback of using excited states, such as 
for quantum information processors, quantum simula¬ 
tors, and quantum metrological applications. 
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Appendix - In order to evaluate (1131) . we will write the 
coupling g k in the electric dipole approximation [H, [Hj 

as g k = -i\J 2 hule 0 A k( r ) ' d where A k (r) defines the 
vector potential in electric dipole approximation, and d is 
the dipole operator for the transition. In terms of vacuum 
Rabi frequency go = 2i^J-^-A k (r) ■ d, we can write the 

coupling as gk = — yjjjju Changing the summation in 

equation (1131) to integration over fc, and using the form 
of gk we write 


m = 


V" 2 o 

167t 3 A 2 


IdO |<7 0 1 4 dk k 2 \a k 


2 1 — COS U) k t 


Integrating over solid angle dO , we write JdO |<7o| 4 = 
where d is the dipole moment of the transition. 
Now changing the variables of integration in the above 
equation to u> using the dispersion relation for free space 
u) = ck , and also using the form of a(u) the equation for 
the decoherence function r(f) simplifies after integrating 
using method of residues: 

_ a4|aol 2 lrf| 4 A r A t 

U 57r 3 e 2 n 2 c 3 A 2 L 2 + A 2 

(wq — A 2 )(l — e~ xt cosa>o t) — 2w 0 Ae _At sinwoU 
+ 2 (wq + A 2 ) 2 J 


Defining the Rabi frequency of the laser as 11 = 
an< ^ the spontaneous decay rate [l5| as T s = 

V ^ 2 ^ 3 0U! ■ Therefore, we can have the following quan¬ 
tity = ^ 4^2 M 4 . Averaging this value over the 

solid angle dO we get = 57 ^ p e ^ 2 e iv A 2 Using this 

expression in the above equation we obtain the decoher¬ 
ence function. Now averaging the function over the 
solid angle and using the result (|go| 2 ) = we get 

the standard result T s = [HI: where we have set 

lo = UJq. 
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